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Abstract: The J + associated production at the LHC is an important process in 
investigating the color-octet mechanism of non-relativistic QCD in describing the processes 
involving heavy quarkonium. We calculate the next-to-leading order (NLO) QCD correc- 
tions to the J/ip-\- associated production at the LHC within the factorization formalism 
of nonrelativistic QCD, and provide the theoretical predictions for the distribution of the 
J/ij) transverse momentum. Our results show that the differential cross section at the 
leading-order is significantly enhanced by the NLO QCD corrections. We conclude that 
the LHC has the potential to verify the color-octet mechanism by measuring the J /ijj + Z^ 
production events. 
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1. Introduction 

Since the discovery of the first charm-anticharm bound state J/V', the study of quarkonium 
is a very interesting topic in both a theoretical and an experimental point of view, which 
provide a good place to probe both perturbative and non-perturbative aspects of QCD 
dynamics. 

In the early days of quarkonium physics, the production and decay of a quarkonium was 
described by the color-singlet mechanism (CSM) [1-5] with the color-singlet QQ pair having 
the appropriate spin, angular-momentum and charge conjugation quantum numbers, and 
it was assumed that the decay and production of a quarkonium can be factorized into 
a short distance part, which can be calculated in a perturbative series of the running 
coupling constant as(M), and a long distance part, which relies on the nonperturbative 
dynamics of the bound state. However, the CSM has encounter many difficulties in various 
theoretical [6-11] and experimental aspects [12,13], such as the appearance of a logarithmic 
infrared divergence in the case of NLO P-wave decays into light hadrons and the huge 
discrepancy of the high-p^ J/il) production between the theoretical prediction and the 
experimental measurement at the Tevatron. 

In 1995, the nonrelativistic QCD (NRQCD) [14], proposed by Bodwin, Braaten and 
Lepage (BBL), provides a rigorous theoretical framework for the description of heavy- 
quarkonium production and decay. In the NRQCD, the idea of perturbative factorization 
is retained, the process of production and decay of heavy quarkonium is separated into 
two parts: short distance part, which allows the intermediate QQ pair with quantum 
numbers different from those of the physical quarkonium state, and the long distance matrix 
elements (LDMEs), which can be extracted from experiments. The relative importance of 
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the LDMEs can be estimated by means of the velocity scaling rules [15]. If one only retains 
the lowest order in v, the description of S-wave quarkonia production or annihilation 
reduces to the CSM. In the case of P-waves, infrared singularities which appear in some 
of the short distance coefficients can be absorbed into the long distance part of color-octet 
S-wave states. Including the contribution of color-octet states, the large discrepancies 
between the experimental data of J/^p production at the Tevatron and the theoretical 
predictions based on CSM successfully reconciled [16-19]. 

Recently, substantial progress has been achieved in the calculation of high order QCD 
corrections to J/tp hadroproduction in order to clarify the validity and limitation of the 
NRQCD formulism. The DELPHI data more favor the NRQCD color-octet mechanism 
(COM) predictions for the 77 J/ip+X process [20,21]. Similarly the recent experimental 
data on the J/i/j photoproduction of HI [22] are fairly well described by the complete NLO 
NRQCD corrections [23], and give strong support to the existence of the COM. However, At 
B-factories, a series of processes was calculated up to the QCD NLO corrections in the CSM 
[24 34] . Together with the relativistic correction [35-39] , it seems that most experimental 
data could be understood. Additionally, the J/tp polarization in hadroproduction at the 
Tevatron [40, 41] and photoproduction at the HERA [42] also conflict with the NRQCD 
predictions. Therefore, the J/ip production mechanism is still a big challenge. The further 
tests for the CSM and COM in NRQCD are still needed on the production of heavy 
quarkonium. 

In order to investigate the effects of the COM in heavy quarkonium production, it is an 
urgent task to study the processes which heavily depend on the production mechanism. The 
J/V" production associated with a gauge boson at the LHC, is a suitable process for studying 
COM. In high energy collider experiments, W^, and J/tp can be identified by using their 
purely leptonic decays [43,44], which are particularly useful in hadron colliders, because 
they provide an enormous suppression of the background. In reference [45], the authors 
give a theoretical prediction at the LO for associated production of heavy quarkonim and 
electroweak bosons at hadron colliders. The numerical results show that for J/ip + ^ 
associated production the CSM provides the main contribution, but for J/ip + W and J/ip 
+ associated production, the COM contribution to the cross section is dominative. The 
NLO QCD corrections to pp J/ip + 7 production in the CSM were provided in [46], 
and the results show that the cross section in large pr region of J/V' can be enhanced by 
two orders in magnitude. Recently, the complete NLO QCD corrections to J/ip + W were 
calculated at the LHC in nonrelativistic quantum chromodynamics [47]. There only the 
color-octet contributes to the process pp ^ J/ip + W at both LO and QCD NLO. The 
numerical results show that the differential cross section at the LO is significantly enhanced 
by the NLO QCD corrections. In this paper we calculate the NLO QCD corrections to the 
associated J/ip production with a Z^ gauge boson in the NRQCD at the LHC, and present 
the theoretical predictions for the px distribution of J/tp. 

This paper is structured as follows: In Sec. II we give the calculation description of 
the LO cross section for the pp — t- J/ip + Z^ + X process, and the calculations of the NLO 
QCD corrections are provided in Sec. III. In Sec. IV we present some numerical results 
and discussion, and finally a short summary is given. 
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2. LO calculations for the — t- J/^ + Z° + X process 



At the LO, there involves two types of partonic processes, which contribute to the pp 
J/ip + + X process: 
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cc[n] + Z^, n 



cc\n\ + Z ,n 



1q(8) 3r,(8) 3p(8) 
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(2.1) 
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where q represents all possible light-quarks(u, d and s). The tree-level Feynman diagrams 
are shown in Fig.|l]. The cross section for the production of a cc pair in a Fock state 
n, a[ij —7- cc[n\ + Z^], is calculated from the amplitudes which are obtained by applying 
certain projectors onto the usual QCD amplitudes for open cc production. In the notations 
of Ref. [48]: 



A 



Tr 



g=0 



Ci/gn?^ 



g=0 



•^CC-[3p(8)] 



g=0 



where ^ denotes the QCD amplitude with amputated charm spinors, the lower index q 
represents the momentum of the heavy quark in the QQ rest frame. Ho/i are spin projectors 
onto the spin singlet and spin triplet states. C^/g are color projectors onto the color singlet 
and color octet states, and £a and Eap represent the polarization vector and tensor of the 
cc states, respectively. 

Then the LO short distance cross section for the partonic process i{pi)j{p2) — ^ cc[n] (^3)+ 
Z^{p4) is obtained by using the following formula: 



a [ij cc[n] + Z°) 



1 



dt^\ALO?', (ij = uu,dd,ss,gg). (2.3) 



The summation is taken over the spins and colors of initial and final states, and the bar 
over the summation denotes averaging over the spins and colors of initial partons. The 
Mandelstam variables are defined as: s = {pi + P2)'^,t = ipi — PzY-,u = {pi — PiY- 
The LO cross section for the pp — )• J/^ + Z^ + X process is expressed as 



a [pp 



J/il) + Z° + X) — 3_ I dxidx2a {ij cc[n] + Z° 

1 + Oij J 



NcoiNpoi 



(2.4) 



where Gij/j^ ^ are the parton distribution functions (PDFs). A and B refer to protons at 
the LHC, Ncoi and Npd refer to the numbers of colors and polarization states separately [48] . 
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Figure 1: The tree- level Feynman diagrams for the partonic processes qq ^ J/tp + Z° (1-4) and 
gg^J/i^ + Z° (5-12). 



The hadronic matrix elements < On^'^ > are related to the hadronization from the states 
cc[n] into J/ip which are fully governed by the nonperturbative QCD effects. The indices 
i,j run over all the partonic species and n denotes the specific color, spin and angular 
momentum cc-state. 

Our numerical calculation shows that the main contributions come from the gg — > 
ccpsj;"'^'*] + and qq — t- cc['^S^^] + Z^ partonic processes. The contributions from other 
partonic processes are relative small. Therefore, we consider only the NLO QCD corrections 
to the gg — > ccpsj^^-*] + Z^ and qq — )• ccpsj^''] + partonic processes in the following NLO 
calculation. 



3. NLO QCD corrections to the pp — > J/t/j + + X process 

In our NLO QCD calculations we consider the following contribution components for the 
pp ^ J/i/j + Z^ + X process: 

► the virtual corrections to the partonic process g{pi)g{p2) cc[^S^^]{ps) + Z^{pi) 
and q{pi)q{p2) ^ cc[^sf\pi) + Z\pi). 

► the real gluon emission partonic processes g{pi)g{p2) ^ cc^'S^i\{p:j,) + Z^{p^+g{p^), 
q{pi)q{p2) -> cc[^s[\p3) + Z^{p^)+g{p^) and q{pi)q{p2) ^ ccfs[\p3) + Z°ip^) + g{p5). 

► the real light- (anti)quark emission partonic processes q{q){pi)g{p2) cc^ s'^'']{p^) + 
Z'^ipA) + g(g)(P5) and q{q){Pi)9{P2) ^ cc[^sf]{ps) + Z\p^) + q{q){pr,). 

► the coUinear counterterms of the PDFs. 



3.1 virtual corrections 

The QCD 0{as) virtual corrections come from the one-loop diagrams including self-energy, 
vertex, box, pentagon and counterterm diagrams. Some representative one-loop Feynman 
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diagrams are displayed in Fig.^. We use FeynArts to generate the diagrams for all related 
subprocesses [49] . Then we use our in-house Mathematica program to simplify and square 
the amplitudes. The phase space integration is implemented by applying FormCalc pro- 
grams [50]. There are UV, IR and Coulomb singularities in the the virtual corrections. We 
adopt the dimensional regularization (DR) scheme to regularize the UV and IR divergences, 
and the modified minimal subtraction MS and on-mass-shell schemes to renormalize the 
strong coupling constant and the quark wave functions, respectively. We adopt the defini- 
tions of one-loop integral functions in Ref [51,52], and use the Passarino-Veltman reduction 
formulas to reduce the tenser integrals to scalar integrals. 




(1) (2) (3) Z 




Figure 2: Some representative QCD one-loop Feynman diagrams for the partonic processes 
gg cc^S^l\ + (1-6) and qq cc[^s[^^] + (7-12). 

We use the expressions in Ref. [53] to deal with the IR singularities in Feynman inte- 
grals, and apply the expressions in Refs. [54-56] to implement the numerical evaluations 
for the finite parts of N-point scale integrals. In the virtual correction calculations, we find 
that Fig.|2|(6,9,ll) contain Coulomb singularities, which are regularized by a small relative 
velocity v between c and c [57], and are canceled by those stemming from the NLO QCD 
correction to the operator < On^^ >■ 

3.2 Real gluon/light-(anti) quark emission corrections 

For the real emission partonic processes, there are IR singularities in the phase space 
integration. We use the two cutoff phase space slicing method (TCPSS) to perform the 
integration over the phase space of these real emission partonic processes [58] . In our calcu- 
lations, the real gluon emission correction contains both soft and collinear IR singularities, 
which are isolated in soft gluon region(£^5 < (5s\/I/2) and hard gluon region(ii^5 > Ss^/I/l) 
respectively. It is easy to find that soft singularities caused by the diagrams emitting 
soft gluon from the charm-quark and anti-charm-quark respectively in J/ijj are canceled 
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Figure 3: Some representative Feynman diagrams for real gluon/light-(anti)quark emission par- 
tonic processes. 

with each other. The hard gluon region is divided into the hard cohinear region (HC),i.e., 
— ti5(or — 125)< ^cS, and the hard noncoUinear region {HC), i.e., — ti5(and — 125)> ^cS 
where ti5 = {pi — p^)^ and t25 = {p2 — Ps)^- Each of the real light- (anti) quark emis- 
sion processes contains only coUinear IR singularity, and can be dealt with in the hard 
collinear region too. In the HC region, the real emission corrections are finite and can be 
calculated numerically with general Monte Carlo method. After summing the virtual and 
real gluon/light- (anti) quark radiation corrections, the remained collinear divergence can 
be cancelled by the collinear counterterms of the PDFs. By adding all the contributions 
together, a finite total cross section is obtained. 

The cross sections for the partonic processes gg — >■ cc[^s[^'']+Z^+g and qq — >■ ccpS'^^^]+ 
+ g m. the soft region can be expressed as 

(gg ^ ccp5«] + + 5) = -g • (-3) • g{p,,P2)aLO?S^^\ (3.1) 



a' {q-q^c-cfsf^] + Z'+g)=-^ 



1 7 



--^{9{Pl,Pc)+9{P2,Pc)) 

where g{pi,Pj) is soft integral function defined as [59-62] 



9{Pi,Pj) = 



(27rM) 



2e 



d^-'p. 



27r Je5<SsVs/2 



Pi 



^jPiPj) 

{piPb){pjPb) {piPbY {pjPbY 



(3.2) 



(3.3) 



while the cross section of partonic process qq — t- ccpS'j;"'^''] + Z*^ + is finite and can be 
evaluated directly in four-dimension using Monte Carlo method. 
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The collinear counterterms of the PDFs, 6G,i/p{x, iif) (P = p ; i = g,u,u,d,d,s,s) 
contains two parts. One is the colhncar gluon emission part SG^^^p°^\x, fXf), another is 

the cohinear hght-(anti)quark emission part dG^J^^''\x, fif). Their analytical expressions 
are presented as follows. 



^Gq{g)/p{x, nf) = SG^^lg^li {x, /X/) + (^Gg^g)"/P {x, fif), {q = u, u, d, d, s, s), (3.4) 
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1 dz 



d,d,s,s 



Pqq{99)i^)Gq{g)/p{x/z,Hf), 



^qg 

{z)Gg/p{x/z,IIf), 



/ -P,,(^)G,/p(x/z,/x/).(3.5) 



q=u,u 



The soft and collinear IR singularities from virtual corrections can be canceled partially 
by adding with the contributions of the real gluon/light-(anti)quark emission processes and 
the gluon emission part of the PDF counterterms ^G^q(^^°^ ■ And the remaining collinear IR 
singularity can be canceled by the contributions of the collinear light-(anti)quark emission 
part of the PDF counterterms 6G^^f^yp exactly. All these cancelations are verified nu- 
merically in our numerical calculations. The explicit expressions for the splitting functions 
(U = qQ,Qg,gQ,99) can be found in Ref. [58]. 

3.3 NLO QCD corrected cross sections 

The NLO QCD corrected hadronic cross section for J/'tp + Z*^ associated production at the 
LHC can be written as 

^QCD ^ ^0 + ^^QCD ^ ^0 + ^^gg 



35(1) 35(1) 



11 + Aa^" 

1 ' 



35(8) 



(3.6) 



The NLO QCD correction Aa^^^ contains three components: (1) Aa^^^i), the NLO 
QCD correction to the subprocess gg ccpS^^^] + (2) Aof the NLO QCD cor- 
rection to the subprocess qq ccfsf^\ + ZO, (3) Act'^^(8), the NLO QCD correction 



to the subprocess 



—7- cc^S^^\ + . Due to the absence of the LO cross section 
for the qq cc[^5|^^] + the A(T^^(g) contains only the contribution of subprocess 



The AcTg^j^) and A(Tg'^(g) can be expressed as: 



35(1) 



{99) _|_ ^{99)9 j_ „{99)<1 _L „{99)<1 
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'35(8) 
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'R 
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'R ' 



(3.7) 
(3.8) 
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where ay''^ , o'^"'^^, c^"'-''^ and f^"'''^ {ij = gg,qq) represent the virtual, the real gluon 
emission, the real light-quark and light-antiquark emission corrections to the cross section, 
respectively. 



4. Numerical results and discussion 

As a check of the correctness of our calculations, we compare our LO numerical results 
with the previous work in Ref. [45]. We adopt the same input parameters as in Ref. [45], 
and reproduce the LO distribution of p'^^^ which is in good agreement with that shown in 
Fig.5(a) of Ref. [45]. 

In the following numerical calculations for the — )• J/^ + + X process at the LHC, 
we take CTEQ6L1 PDFs with the one-loop running as in the LO calculations and CTEQ6M 
PDFs with the two-loop as in the NLO calculations [63]. For the QCD parameters we take 

(3) (3) 

the number of active flavor as rij = 3, and input Aq^j-, = 249 MeV for the LO and Aq^j-, = 
389 MeV for the NLO calculations [23], respectively. The renormalization, factorization, 
and NRQCD scales are taken as fir = fJ-f = rriT and = "ic, respectively, where rriT = 

(pt^^ + "^j/^ is the J/V' transverse mass. The masses of the external particles and 
the fine structure constant are taken as rriz = 91.1876 GeV, rric = mj/^/2 = 1.5 GeV and 
a = 1/137.036. The ( C'"^/V'[3pj8)j ^ ^ ^^-^^2) LDMEs satisfy the multiplicity relations 

( O-^/^i^pf] ) = (2J + 1)( 0^/'^[3P(S^)] ), 
and the linear combination exists as 

We take the LDMEs ( ©■^/'^[^^f^] ) = 1.3 GeV^, ( 0^/^[^s[^''] ) = 2.73 x IQ-^ GeV^, 
Mr = 5.72 X 10^ GeV and r = 3.54 as the input parameters, which were fitted 
to the Tevatron RUN-I data by using the CTEQ4 PDFs and taking into account the 
dominant higher-order effects due to the multiple-gluon radiation in the inclusive J/ip 
hadroproduction [64]. Then we fix the ( ©■^/'^[^sj^^] ) and ( O-^/'f'l^PP] ) LDMEs by the 
democratic choice [65, 66] 



2 

We have checked numerically the independence of the total NLO QCD correction on 
the cutoffs 6s and 6c in the range 1 x 10~^ < (5s < 1 x 10~^ and 6c = (5s/50. The independence 
exists in our numerical results within the error tolerance. In further calculations, the two 
phase space cutoffs are fixed as 5^ = 1 x 10^^ and 6c = 6s/50. Considering the validity of 
the NRQCD and perturbation method, we restrict our results to the domain p^^^ > 3 GeV 
and \yj/^\ < 3. 

The dependence of the cross section on the renormalization scale fj,r and factorization 
scale fif induces uncertainty for theoretical prediction. In Fig.^ the /i dependence of the 
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Figure 4: The dependence of the LO and the NLO QCD corrected cross sections for the process 
pp — >■ J /il> + Z'^ + X on the factorization scale and rcnormalization scale (/i//io) at the LHC where 
we define n — — Hr and = rriT ■ 



LO and the NLO QCD corrected cross sections with the constraints of p'l^^ > 3GeV and 
\yj/i>\ < 3 for pp —7- J/'ip + Z^+X process, is shown with our default choice fi = Hr = fJ-f and 
the definition of /io = nT-T- There fx varies from /io/3 to S/Uq. Not like the usual expectation, 
Fig.^ shows that the NLO QCD correction can not improve the LO scale independence for 
the pp ^ J/ip + -\- X process at the LHC. Actually, the similar behavior appears also 
in the results at the NLO in NRQCD (see Ref. [19]). The related theoretical uncertainty 



amounts to at the LO and to _2i% at the NLO when goes from /io/2 to 2fj,Q. 

In Fig.||, we present the LO and NLO QCD corrected distributions of p^^^ and the 
corresponding K-factors for the pp — t- J/tp + Z^ + X process at the LHC. For comparison, 



we also depict the LO and NLO QCD contributions from the ccpsj and cci'^S'i'] Fock 
states and their corresponding K-factors in these figures separately. The upper figure shows 
that the differential cross sections at the LO are significantly enhanced by the NLO QCD 
corrections and both the LO and NLO QCD corrected differential cross section curves drop 
rapidly with the increment olp'^^ . We can read out from the figure that the LO differential 
cross section of p'Jj!'^ decreases from 79.2 fb/GeV to 0.3 fb/GeV and the NLO QCD 
corrected differential cross section oi p'^^ decreases from 221.1 fb/GeV to 1 fb/GeV as 
p^^^ goes up from 3 to 50 GeV. The corresponding iT-factor, defined as K = / '^"j'j'^ , 

varies in the range of [2.79, 3.42], and reach its maximum when p"!^"^ = 50 GeV . 

From Fig.^, we find that the contributions to the differential cross section from the 
cc[^5^^^] and cc[^S'^'^^] Fock states are comparable at lower p'^J'^ . As the p'^^ increases, the 



:r3c(8)i 



contribution of cc[^5|^^] decreases much faster than cc[^5j'''']. In all the plotted p'^rj!^ range, 
the CO contributions to the differential cross section of the process of pp — >• J /ip+Z^+X are 



■3 0(8)1 
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Figure 5: The LO and NLO QCD corrected distributions oip'^J^ and the corresponding K-factor 
for the pp —T' J/ip + + X process at the LHC. 

dominant, and that will be beneficial to our study of CO mechanism. For the contribution 
from the pp — )• cc process, the ii'- factor of differential cross section increases 

from 0.55 to 1.50 as p'!^'^ increases from 3 GeV to 50 GeV. As for the pp — )■ cc[^5j^^]+Z''+X 
process, the ii'-factor is larger than that for the pp — t- ccpS^"*^^] + Z*^ + X process, which 
decreases from 4.9 to 3.6 as p'^^^ increases from 3 GeV to 50 GeV. 



5. Summary 

The J/tp + Z^ associated production is an important process in investigating the production 
mechanism of J/V'- In this paper we investigate the NLO QCD corrections to the J/ip + Z^ 
production at the LHC, and present the numerical predictions of the p'^^^ distribution of 
the p^^^ up to the QCD NLO. We find that the LO differential cross section for the process 
pp ^ J/tp + Z^ +X is heavily enhanced by the NLO QCD corrections, and the K-factor can 
reach the value of 3.42 in the large p'^^^ region. The NLO QCD corrected differential cross 
section can reach 221.1 fb/GeV in the vicinity of p^^'^ ~ 3 GeV. The process of J/ip + Z 
associated production is dominantly contributed by the CO mechanism in the large p^^^ 
range. We conclude that the LHC has the potential to study the J/ip + Z^ production 
process, by which we can investigate the production mechanism of heavy quarkonium and 
extract the universal NRQCD matrix elements. 
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